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UDC 593.3 


STUDY OF DYNAMICS OF CIRCULAR RING STRUCTURES ON PROGRESS-40 CARGO 
SPACECRAFT 


917F0158A Kiev DOKLADY AKADEMII NAUK UKRAINShOY SSR: SERIYA A -- 
FIZIKO-MATEMATICHESKIYE | TEKHNICHESKIYE NAUKI in Russian No 11, 1990 
(manuscript received 14 Aug 90) pp 36-43 


[Article by academician B. Ye. Paton, Yu. P. Semenov, P. M. Belousov et 
al., Kiev Construction Engineering Ins*°itute] 


[Text] A space experiment to unrol! two large structures 20 m in 
diameter each, which use drives from material having a memory effect for 
unrolling, was conducted on the Progress-40 spacecraft and the Mir 
orbital station during the period from 3 through 5 March 1989. 


The process of unro!ling and shaping these structures under space f] ight 
conditions was studied experimentally, and the stiffness and dynamic 
characteristics of the unrolled loop structures during dynamic 
operations with the cargo spacecraft were studied. 


The free vibrations of large circular structures and their transient 
motions, caused by programmed progressive and rotational maneuvers of 
the carg» spacecraft, were calculated in preparation of the experiment. 


Let us consicer the problems of mathematical modeling of the dynamics of 
the carrier--ing system, the center of mass of which moves in a circular 
orbit about the earth. When performing maneuvers of the carrier 
according to the given program, the rings connected to it perform 
complex motion, consisting of translational motion of the body itself 
and the relative motion of the rings with respect to it. Let us 
disregard the orbital motion caused by gravity and inertia and let us 
assume that the system is free to study the dynamics of relative motion 
of the rings. Let us introduce the inertial coordinate system 0,XY7 and 
the system O.xyz, rigidly bound to the carrier so that its axes coincide 
with the principal central axes of the system, and axis O.7 passes 
through the centers of rod rings, while axis O-x is perpendicular to 
their plane. 


The equations of dynamic equilibrium of the deformed element of the rod 
ring are represented in the form [4] 











; dN 
Few F-t-0, 


73 7. w<M--mF= 0 (1) 
where F is the internal force vector, M is the internal moment vector, s 
is a coordinate measured along the axis of the line of the rod by the 
distance from some initial point to the current point, f is the vector 
of the external load intensicy, which also includes the force of 
inertia, fT is the unit vector of the tangent to the axia! line of the 


rod, # is a Darby vector, and d/ds denotes the local derivative. 


Let us introduce the coordinate system (u, v, wW), bound to the 
considered cross-section of the rod, axis w of which is directed along 
the tangent to the flexible ine, and axes u and \ are directed along 
the principal central axes cof inertia of the cross-sectional area. 


Projections My, My, and M, of moment M the corresponding axes are 
expressed by the stiffnesses upon bend f A and B and by twisting of 
C, and also of curvature p, q of elem s and twisting r of the rod 


axis by the relations 


in which the values of p, q, andr a deformed state of the rod are 
dentoed by zero. 


Let us supplement system of equatio..s (1) and (2) kinematically by 
Frenier equations [1] 
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(3) 





Here n and b are unit vectors of the normal and binormal and p and T are 
the radii of curvature and twisting. 


Let us also write the equations 
dx dy 


=m = . (4) 


—— \ . 
ds x as us 


The system of solving equations (1)-(4), formulated in the domain of 
variation 0 < s ¢ S of independent variable s has a common 18th order. 








The first six integrals, which follow from the conditions of 
orthonormalization of the Frenier basis, 


close the system of solving equations and permit us to avoid 
redetermination of it. 


To construct the vibration equations of the rod, let us use the 
D’Alambert principle and let us represent the force acting on it in the 
form 


f - va, (5) 


where 7 is the linear density of the rod and a is the acceleration 
vector of the element of the rod in system O,XYZ. 


Let us use the Coriolis formula [2] 4=a.—a,— a, where a, = xi + yj ok 


is relative acceleration, a =a,—+ex.p—2»(Qxp) is translational 
acceleration, a 20> (xi — yj — zk) 18 Coriolis acceleration, ap is the 


acceleration of the center of mass of the carrier, f and ¢€ are vectors 
of the angular velocity and acceleration of the carrier, i, j, and k are 
unit vectors of system O.XYZ, and p = xi + yj + zk; differentiation with 
respect to t is denoted by the dot. 


Since vectors a,, a,, and a., used for calculation of f, are calculated 
in coordinate system O,XYZ, we initially find its projections 


from which by formulas 


we turn to components f,, fy, and fy, used in equations (1). 








Having substituted ordinary derivatives for independent variables s and 
t in equations (1)-(4) by partial derivatives, we find the solving 
system of nonlinear equations of relative vibrations of the rings. 


When studying the dynamics of the entire system, one must construct the 
equations of motion of the carrier and solve them jointly with the 
equations of motion of the rings. To do this, let us use the theorems 
on variation of the momentum Q and of the moment of momentum K of the 
carrier in the forms 


Q R - \i vi (6) 


where R and M* are the principal vector and principal moment of 
reactions at the points of finishing of flexible rings, caused by their 
elastic deformation upon vibrations of the system, and F? and M® are the 
principal vector and principal moment of all active forces. 


Using equations (1)-(6), the dynamics of the induced and free motions of 
the system is studied. Let us dwel! on the results of solving the 
problems of free oscillations. 


Let us limit ourselves to calculation of the lower values of the range 
of vibration frequencies and to construction of the corresponding forms. 
This permits one to turn from the initial continual problem to the 
corresponding discrete problem and to consider only a finite number of 
degrees of the system. 


Let us assume that the vibrations of the system are determined by 
time-dependent displacements of the system of its discrete points 
where N is the number of discrete points. The 


equations of free motions of the system can then be represented in 
matrix form 


WX -—AX- 0 (7) 


where X is the vector of displacements of the nodal points of the 
system, M is the matrix of masses, and K is the matrix of the stiffness 
of the discrete elements. 


Matrix M is determined by the mass of each extended element, 
concentrated at the nodal points. 


The stiffness matrix K can be determined by using the property K = $°!, 
where ® is a pliancy matrix. The eleme:its of matrix $ are the 
displacements alon’ the permissible couplings of the directions upon 








action of the corresponding uiit concentrated force on the node [6]. 
They are found as a result of solving the problems of the dynamic 

equilibrium of the rings upon corresponding steady translational and 
rotational motions of the system [4]. 


\ 


Assuming that free motion of the system has the form of won re 

vibrations, let us represent the solution of the system in the form 
X(t) = H sin wt, where H is tne vector of the form of vibratioré and 
is the circular frequency corresponding to it. 
solution into (7), we find 


The eigenvalues @#, and the eigenvectors Hj 


QR-algorithm. 


(K — wo-M) H — 0. 


Substituting this 


(8) 


are calculated using the 








J) te GO toe 


1. Number of form 


2. Period, s(Tj ) 
3. Frequency, rad/s (@;) 


0357 
0474 
, 1006 
,1102 
1566 
2340 
» 2048 
, 4016 
4398 
4408 
4411 
»4439 


Angle of rotation, rad 


Displacement in direction of 


axis 


A total of 18 values of the spectrum of the natural vibration 

frequencies of the carrier-elastic ring system is calculated by the 
developed method, and the corresponding forms of motion are constructed. 
The results of calculation are presented in the table. 


The ordinal 








numbers of the eigenvalues are given in the first column, the length of 
the vibration period T; is given in the second column, the corresponding 
values of the circular frequencies «, are given in the third column, the 
amplitude values of the angles of rotation of the carrier with respect 
to the corresponding axis are given in the fourth column, and the axis, 
in the direction of which,the greatest displacement of the typical point 
of the ring is realized, ys indicated in the fifth column. The forms of 
Vibrations of the rings are normalized so that the amplitude value of 
the maximum displacement is equal to one. 


The axonometric images of the first six forms of natural vibrations are 
shown in Figure 1 with respect to a fixed coordinate system. Two 
coordinate systems--a fixed system OXYZ and coordinate system Oxyz, 
rigidly bound to the carrier, rotated during its vibrations with respect 
to the fixed system--are also presented in it. The movable coordinate 
system does not happen in cases when the rotational motion of the body 
is absent. The forms of vibrations of the rings are normalized in all 
figures so that the value of the maximum relative displacement is equal 
to 1/5 of the diameter of the ring. 


As calculations showed, the vibration modes that permit rotation of the 
carrier have higher frequencies, compared to the descending modes of 
natural vibrations of the rings, rigidly fastened to a fixed support. 


Upon study of the dynamic and strength characteristics and also of the 
controllability of the cargo spacecraft with unfolded large annular 
structures, the following problems were solved: 


experimental determination of the frequency and decrement of the main 
tone of vibrations of the structures; 


2) a check of the capability of the structures of restoring a given form 
after the cargo spacecraft completed turns and correction pulses; 


a check of the results of preliminary calculations of the structural 
dynamics and applicability of the calculating models to be used. 
The most complete information, including a television image of the 
vibrations of one of the annular structures, transmitted from the cargo 


spacecraft igure 2), during the experiment was found in two sessions 
n which the cargo spacecraft accordingly performed a correcting pulse 
AV = 0.6 m/s and a programmed yawing turn by angle » = 120° at angular 
velocity «# = 0.67 /s. 


Symmetrical vibrations of the structures were excited with perioc 190 s 
during the correcting pulse. The points of the rings most remote from 
the carrier were deflected from the neutral position by 11 m. These 
values of the period and amplitude of vibrations essentially coincided 
with the results of earlier calculations. It should be noted that the 
actual period of vibrations was determined immediately from three 








sources: by a television image from the cargo spacecraft, oy a report 
of the crew of the "Mir" spacecraft and by telemetry data. 





Figure 1 


The dissipative properties of the large annular structures were studied 
upon analysis of the vibrations caused by the correcting pulse, which 
are difficult to determine by calculation or by ground tests. Observing 
these vibrations directly, the crew of the "Mir" orbital station noted 
Slow damping of them, while preliminary analysis of the telemetry data 
permitted an estimation of the logarithmic decrement of 4 = 0.05. 








Figure 2 





When considering the methods of improving the ntrollability of a 
spacecraft with large unfolded annular structures and taking int 








account the weak damping of vibration of these structures, main 
attention should be devoted to dynamic suppression of vibrations, 
excited during maneuvers of the spacecraft, rather than to passive 
damping. Thus, the dynamic effect of essentially total instantaneous 
correction of vibrations of the structure as a result of selecting the 
moment of stopping the rotacion of the spacecraft can be specifically 
used as one of the principles of selecting its optimal control. 





Figure 4 


The dynamic effect of essentially total kinematic suppression of 
vibrations was demonstrated when the Progress-40 spacecraft performed a 
programmed turn, and the angle of rotation of 9» = 120° was previously 
calculated on the basis of the measured value of the actual period of 
the main tone of vibrations of the annular structures--T = 190 s and of 
the nominal rotational velocity of the cargo spacecraft w = 0.67°/s, 
with the calculation that the structure performed a single vibration 
during the turn of the cargo spacecraft. 


It was obvious from the television image from the Progress-40 spacecraft 
(Figures 2-4) that the annular structure performed one vibration during 
the turn, and that it essentially stopped in the neutral position at the 
moment of completion of the turn and the beginning of stabilization of 
the cargo spacecraft. The theoretical conclusion of the possibility of 
rapid suppression of vibrations, having the value of the optimal control 
of the dynamics of large structures when working out the methods, was 
thus confirmed. 








Comparison of the results of calculation and data of full-scale tests 
indicates that they are in satisfactory agreement. 
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[Article by N. V. Vasilenko, S. A. Sarapulov, and A. M. Pavlovskiy, Kiev 
Polytechnical Institute] 


[Text] The effect of progressive and angular vibration on the motion of 
traditional mechanical gyroscopes has now been studied. However, the 
problem of the effect of vibration of the base on the accuracy of the 
new inertial information meter--the solid-state wave gyroscope (TVG), 
which is of considerable interest, has been insufficiently studied 
[1-4]. 


The nonlinear dynamics of a solid-state wave gyroscope with imperfect 
Cavity on a progressively vibrating base is considered in this paper. 


The state of the object is given by generalized coordinates pp(t) and 
Gn(t), corresponding to conjugate natural forms of vibrations of a 
hemispherical cavity Unn(8)cosnq, Urn(9)sinn@ (kR=1, 2, 3),, which are 


introduced by the approximation 


te ’ 
— . 


n 


where angle » is counted in the circumferential direction and angle @ is 
formed by the axis of symmetry and by the internal normal. Application 
of the Bubnov-Galerkin procedure to the nonlinear equations of motion of 
the cavity on a vibrating base [4] with regard to representation of 








defects of a hemisphe. ical dome in the form of expansions into a Fourier 
series by circumferential angle » [5] results in ordinary differential 
equations in generalized coordinates that describe nonlinear vibrations 
of an imperfect cavity on a progressively vibrating base. Let us 
consider resonant vibrations when the external perturbation frequency y 
is a multiple of the working frequency w2 (the case of excitation of a 
second pair of natural forms of vibrations in the cavity, n = 2, is of 
interest for applied problems [1, 6]), assuming that there is no 
internal resonance in the system. Equations that describe the dynamics 
of the main form have the form 


' _*9 : ‘* ' An Y ’ ol ' L a LL . 7 ate ( aDo 
(Mp + €My,) Jo + EMyoPy + Dek Qn, — hgGg SH (Cg H EC qd HH CyyPe 
+ © (Hoh, = A191) — &f, — Oe) = 9, (1) 
(M, — EMao) Po aa EITM o14o — Qe A ,82G. alana Chg ls im tc, — e¢ 99) J 9 . ( oy4o a 
— &(A pith + Ha - Harty @4 + Oiler) = Q, 


v0 ! 


where 


f, = mWox + mW oy + ni.Woz. 7s nM ? Mt yy —~ mv > 


Here mo, co, Ko and do are the coefficients of inertia and stiffness, 
and gyroscopic and damping coefficients of an ideal cavity, 
respectively, {} is the angular rotational velocity of the cavity about 
its own axis of symmetry (it is considered a slowly variable time 


function, i.e., we disregard terms with {) and (2), m,; and cj; (i, j= 
1, 2) are perturbations of the coefficients of inertia and stiffness due 
to errors in manufacture of the hemispherical cavity (the fourth 
harmonics of expansion of the 


defects into a Fourier series [5]), and 27). ¢).=C2;, My = —My», 
C1;= —Coo: 1, m,* (l=!], 2, 3) are perturbations of the inertial 


characteristics of the cavity, m3 and m3* correspond to the second 
component of the defect in Fourier expansions, m;, m;*, mo, and mo* are 
the first and third harmonic, Wox = —y NS cosvt, Woy = —77 Yo cos 7! 

Woz= —y°Zo cos 1 are projections of the vector of linear accelerations 
of the base onto axes X, Y, and Z of a coordinate system rigidly bound 
to pole O of the cavity, and axis OZ coincides with the axis of symmetry 
of the cavity, q;, pj, and po are generalized coordinates that 
characterize the motion of a hemispherical cavity as a solid 


_ OX COS ¥!. Py — gi) , COS F a7 — bi pyle, \ ) ¥!, 0, CONS, 
: e ! 


‘al 
‘7] 


© vata ? 
£,= const, £, = const, 











where Xo, Yo, Zo and y are the amplitudes and frequency of external 
perturbation, Ho, H;, Hg, and Hg are constant coefficients with 
nonlinear terms, and ¢€ is a small parameter. 


Let us consider two important cases of vibration of the base. 


Longitudinal vibration. The base moves along the axis of symmetry of 
the cavity (Z = Zo cos yt). Analysis of equations (1) shows that the 
following resonance situations are possible in the system: a) y = # + 
+ €9, where ¢ is a detuning parameter. In this case, longitudinal] 
vibration of the base results in excitation of "spurious" induced 
vibrations in the cavity. In zero approximation of the perturbation 
method [7], solution of the system of equations (1) has the form 


gy” (Q) ( T)) exp ( 'W7°)) i Q (T,) exp (--- iWs/'y), 


pi) <= P(T, explo. 4) + P(T,)exp(— iw,7,), 


) 


where 7,=c'f; Q,Q and P, P are complex conjugate amplitudes. Let us 


represent the complex amplitudes in exponential form 


Q = a(T,) exp {i[a(T,) + o7,}}, 
P = )(T,)exp {i (fh (7,) — o7,j}. 


Substitution of variables by formulas 6 + ia = Cexp(2i¥),p—«= 6, 2+/= 


= ¥ permits us to write the condition of the absence of secular terms 
in the solutions of the first approximation q‘!)» and p‘!)» in new 
variables C, ¥, y, and 7. A real instrument has a control system that 
maintains vibrations of the cavity in the form of a standing wave, 
since, as shown in [8], it assumes fulfillment of the following 

Any 


), << - 9. Taking the foregoing into account, 


conditions: ) = 0,—.— 
Jl 1 oT. 


we find a formula that determines the precession rate of the standing 
wave during functioning of a solid-state gyroscope due to the effect of 
longitudinal vibration 


KO 7-Z,, sin (47/2) 


Pe opt + AS (m, cos 2'¥ — m? sin 2) ; 
ym din wl “3 cos ~ 4 Mm, sin ~ P). (2) 
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a(..) 
oT, 


of orientation of the wave pattern, respectively, and ,=(c,m,) is 





Here (...) = C= (b} —.*)'*, Y are the vibration amplitude and angle 


the natural frequency of an ideal cavity. It follows from formula (2) 
that the integrating effect of a solid-state wave gyroscope will occur 
only at angular velocities {}, which exceed the threshold value 


¢) ) P be } 
\mZ w/2K,C. Here Am = (in? — m3°)'?, ¥ W,. 


t) y = 2¢2 + 2€¢. Longitudinal vibration excites “spurious” parametric 
ibrations. The equation that determines angle ¥ has the form 


W’ = —K,Q2m,. (3) 


It follows from equation (3) that longitudinal vibration in the given 
case does not result in an additional "spurious" component in the wave 
pattern. 

Transverse vibration. The base moves in a plane perpendicular to the 
axis of symmetry of the cavity (.\= 15 cos yf, Y=Yovos yt). The following 


resonant si*uations are possible in the system: 


oO 
=™ 
=I 

ho] — 


In both cases, transverse vibration of the base results in excitation of 
“spurious” induced vibrations of the cavity. Expressions that determine 
the position of the wave pattern as a function of the angular rotational] 
velocity and vibration parameters of the base have the form 


ay VW = —(2m,)7' KO = (4m 
=M) KO - (4m w,C) (my X, + md.) cos 24 
— (mA, - M24) sin 24°] sin (4/2), 
1 Ww’ In, \~! . 
db) Y — (2m) AO + (4mwo0y7 | | Hea X45 yr 
09 1k) 4 US 2'F - 
—_ (HM #2 y2 ) 02) 
(Hp) X5 Ayge°Y',) sin QV] sin (9/2) 


In the given case, we will also have some threshold value of the angular 
rotational velocity of the base, below which there is no integrating 
effect. 


It follows from the studies that there are several mechanisms of errors 
in the operation of the instrument upon vibration of the base: 








excitation of additional “spurious” vibrations of the cavity that lead 
to distortion of the readings of the instrument due to the presence of 
errors in manufacture of the cavity; 


excitation of additional “spurious” vibrations due to the presence of 
nonlinear coupling (between generalized coordinates) that characterize 
intermodal interaction. 


The results of the studies can be used in determining the dimensions of 
the zones of absence of the integrating effect of a solid-state wave 
gyroscope, the nominal value of the vibration amplitude of the cavity, 
and also when working out the requirements on the precision of 
manufacture of the cavity and on the vibration protection system. 
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[Text] The problems of increasing the efficiency of utilizing fuel and 
energy resources are acquiring ever more timely value in ali sectors of 
the national economy. Problems of energy conservation will determine tx 
a great extent in the future the development of industry, transport, 
agriculture and the communal-everyday sector. The problem has special! 
significance for Latvia under conditions of transition to independent 
forms of management and republic cost-accounting, since 85 percent of 
the fuel and 50 percent of the electric power are delivered from other 
regions. Along with this, the problem of adequate and real nature- 
conserving activity is valid. Interaction of facilities of the fuel and 
energy complex with the environment is an exceptionally important 
problem. On the one hand, the value of an uninterrupted effective 
energy supply of all sectors of the national economy is very high, and 
on the other hand, the fuel and energy complex has the most significant 
effect on the human habitat. Only centralized power enterprises 
discharge more than 40,000 tons of noxious materials annually in the 
atmosphere in Latvia, which comprises 21 percent of the total discharge 
in the republic. 


The increasing pollution of the atmosphere is assuming a global nature 
today. Discharges of noxious and harmful pollutants into the atmosphere 
by many sectors of industry have a deleterious effect on living 
organisms, destroy the ozone layer, and contribute to the occurrence of 
acid rain. The reality of the greenhouse effect occurring is a 
consequence of variation of the structure of radiant heat transfer as a 
result of intensive utilization of thermal power systems. These 
phenomena may affect the entire country and even whole continents. 
\lthough scientists have different opinions on the time intervals and 
intensity of the ecological changes, nevertheless they are still 
Significant, which requires urgent study. 


4 session of the General Assembly of the Academy of Sciences. devoted to 
discussion of the project of the main propositions of the Program of 
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Biospheric and Ecological Research of the USSR Academy of Sciences ior 
the period up to 2015 in December 1988 noted that we are in the phase of 
anthropogenic effects on nature, when the processes of ecological 
degradation are even more local in nature, but global trends, hazarcous 
for the biosphere as a whole, are already being detected on this 
background. This is the moment when the process of irreversible 
ecological changes can be stopped now through the efforts of people. 


Outside of health and well-being, all our social goals and political and 
ecological ideals are impractical. Alternative versions of solving 

soc ioecological and technical problems must be worked out. Academician 
P. L. Kapitsa convincingly demonstrated the interdependence of the 
growth of the gross national product and energy production for all 
countries of the world. This means that if there is an increase of 
energy, there will also be an increase of the national product. 


The time has come to combine the means of survival into a single 
complex, optimal in the given ecological situation [1]. This is a most 
important scientific task both on a global scale and for each country. 


The Physics and Power Institute of the Latvian Academy of Sciences is 
participating in development of the section "Power engineering” of the 
general academic program of biospheric and ecological research for the 
period up to 2015. The topic being studied is a method of estimating 
the loss due to the effects of energy facilities. Problems of 
estimating the ecological impact of energy-conserving measures, which 
include analysis of the engineering solutions of energy conservation in 
industry and ecological zoning and dipersion of energy resources, 
structural shifts in industry, and their energy and ecological 
consequences, are being considered in the first phase of these studies. 
The studies are directed toward discovering the characteristic features 
of energy conservation with different layouts of energy conservation of 
enterprises, a rational structure of industry and with specific 
principles of estimating the ecological impact of energy-conserving 
plants, now used in Latvia. 


Determination of the principles of ecological estimation of the 
engineering solutions of energy conservation in our republic is based on 
analysis of the characteristic features of natural conditions, 
functioning of the natural and technical environment, and of the state 
of the ecoosystems of the given territory. Thus, optimization 
calculations are made with regard to ecological factors under specific 
conditions of a corresponding region, but the principal approach to the 
problem is identical for all economic regions. 


\nalysis of the functioning of the territorial economic-ecological 
system of the Latvian Republic shows that its intensity determines the 
onsequences of the ecological concept, implemented during the previous 
vears, which was constructed on the basis of hypertrophic development 
f industry, oriented toward extensive growth factors, and was not 
vocompanied by a search for effective ways of reducing the energy and 








resource capacity, and solution of environment protection problems. The 
main industrial centers at present--cities and urban agglomerations-- 
have become areas of extensive negative transformation of the 
environment, the effect of which is very significant not only on the 
nature of adjacent territories, but also on the vast regions of river 
basins, and coastal] waters of the sea. One must sometimes talk about 
the symptoms of a crisis state, when the parameters of the environment 
are approaching the permissible limits of changes, and suppression of 
reproductive capacities and exhaustion of the ecological capacity of the 
territory have been noted; therefore, intensification of the 
environmental protection trend now conducted outside economic reform is 
necessary, the ecology aspect of cost accounting is required, and this 
must be transformed to a reliable barrier on the path toward wasteful 
utilization of nature. 


The current phase of development of the economy is characterized by a 
stable increase of the expenditures of social labor, embedded in 
environmental protection and in reproduction of natural resources. A 
growth of the scales of ecological investment increases the 
responsibility of each planning solution in this field, and imposes 
higher requirements on their effectiveness. At the same time, society 
needs optimal development of TEK [not further identified] as a necessary 
condition of enhancing the life level and national income. An objective 
need arises in this regard to establish interrelationships between the 
results of economic activity and its ecological consequences. 
Substantiation of the economy of production is an inseparable part of 
the management system that affects the selection of priorities in 
support of the national economy with natural and energy resources. The 
main causes of a negative effect of the national economy on the 
environment are imperfection of technology, an increase of the scales of 
production, and its concentration on a limited territory. The increase 
of energy consumption can be covered by 30-40 percent in the Latvian 
republic by introduction of energy-conservation technologies and 
scientifically substantiated standards of energy consumption [2]. 


Life persistently requires extensive analysis of the energy consumption 
of manufactured products. High losses of energy resources are formed 
due to the fact that we do much that is unnecessary, and even in 
weighted form. The unsubstantiated high materials consumption is 
accompanied by impermissibly high energy consumption. Recovery of 
secondary resources that vields significant energy conservation has not 
become widespread. The high degree of depletion of production funds in 
industry, which comprises more than 47 percent, also results in rational] 
consumption of energy in our republic. This part of the morally and 
physically obsolete production funds is also a significant source of 
environmental pollution. 


The fact that tenfold more funds are expended in the country on routine 
maintenance of operated equipment than is required for replacement of it 
indicates the redundancy of this practice [3]. 








Utilization of secondary energy resources can reach more than 46 percent 
in the future structure of measures on conservation of energy resources, 
while the use of energy-conserving technology in industry and in 
agriculture may reach approximately 38 percent [2]. 


An experimental lot of facade heat regulators is functioning 
successfully in Riga. The instruments, developed at FEI, contribute to 
conservation of energy resources and prevent overconsumption of thermal 
energy. The heat loss factor of new multistory buildings is 0.42-0.47 
kilocalories per square meter, which is 25-30 percent higher than the 
optimal value. 


Taking into account the ecological impact is of considerable 
significance in economic estimation of energy-conserving measures and in 
determination of the priority of implementing them with limited monetary 
and material funds. Analysis of the ecological impact is related to 
estimation of the consequences of discharges of energy sources that 
supply energy-conserving plants compared to traditional plants, which 
utilize a closed form of fuel (in the presence of scrubbing structures 
that provide the PDV [maximum permissible discharge} at the required PDK 
[maximum permissible concentration]). The consumption of energy 
resources on the energy supply of the drive of energy-conserving plants 
for use of secondary thermal resources is insignificant, and sometimes 
(utilization of TNU heat pumps) may reach 30 percent or more of the 
extracted alternative energy resources. Therefore, the ecological 
impact is determined by the difference in the amount of fuel burned for 
comparable variants (the drive of plants that conserve energy or 
generate the same amount of energy by the traditional method) and by the 
loss due to discharges. Estimation of the loss is determination of the 
expenditures, related to usefully generated (conserved) energy at 
comparable energy plants (per ton of conventional fuel). The economic 
impact of energy-conserving measures can generally be formulated with 
regard to the ecological loss due to discharges as 


A3’x¢ = Aa + AB, + A37+ Ay, (1) 


where A3 ‘5. is the specific economic impact of energy-conserving 


measures with regard to the ecological impact (rubles/t of conventional 


fuel), 13, is the difference of the specific expenditures, dependent on 


capital investments on a traditional power plant operating on the close 
fuel cycle and on an energy-conserving plant, ‘3. is the difference of 


the specific expenditures, dependent on operating expenditures, 43. is 


the difference of the specific expenditures for fuel and energy for 
comparable variants, and Ay is the difference in the loss due to 
discharges. 














The difference of specific fuel and energy expenditures characterizes 
the conservation of energy resources and can be determined as 





(2) 


l l 
N3y = 25 ( ‘a. a =} (3) 


where Bot and Bsc are fuel consumption for a traditional energy plant 


that utilizes a complete type of fuel and an energy-conserving plant 
(tons of conventional fuel), Qs, and Q5¢ are the useful energy 


consumption (tons of complete conventional fuel), and lo¢ are the 


"37 
energy utilization factors (9,5 is from 0.8 to 0.5 for traditional heat 


supply plants, while fee * 1.0 to 2.0 or more for energy-conserv ing 


plants, for example, for TNU, due to the use of ambient thermal energy), 


and 2,7 are the specific expenditures for a complete type of fuel 


(rubles/ton of conventional fuel). 


The difference in damage due to discharges Ay (upon observation of PDV 
and PDK) is dependent on the difference in the amount of burned fuel per 
unit of useful energy and of the specific damage due to discharges in 
the compared variants, i.e., 





(Bx _ Bx ) (1-1), 
Ay= ~o_ } ¥8= © he (4) 


where Ay is the difference in damage due to discharges (rubles/t of 
complete -onventional fuel) and YR is the specific damage due to 


discharges per unit of burned complete type of fuel (rubles/t of 
conventional fuel). 


When using different types of fuel to drive an energy-conserving plant 
and in a traditional plant, Ay is determined as 








where y are the specific damage due to discharges from 
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combustion of complete type of fuel and of fuel utilized to generate 
energy to drive an energy-conserving plant (rubles/t of conventional 
fuel), respectively. 


Selecting the type of complete fuel and energy to drive energy- 
conserving plants is related to optimization of the fuel-energy balance 
of the region, the structure of which can be changed significantly upon 
implementation of wide-scale measures on conservation of energy 
resources (gassification of a region). To determine the volume of 
conservation of energy resources and of the amount of complete type of 
fuel, one must also resort to using optimization models. 


A method that determines the damage due to discharges Ay is rather 
complex and requires further studies. However, this is a separate 
direction of the studies, related to participation of scientists of 
different specialties on problems of ecology. 


Taking into account the additional economic effect of energy-conserving 
measures in a region significantly affects both their ranking according 
to ecological impact, and the volume of feasible implementation of 
energy-conserving measures for the future. 
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[Text] General Characteristic 


The program system (PK) is designed to study the main directions and 
principles of development of energy associations (EO). The most typical 
problems to be solved by using the ORS E program system are: 1) 
selecting the variants of development of systems-forming networks of 
energy associations at existing voltage by introducing new electric 
power transmission lines or by renovation of them; 2) substantiation of 
the introduction of electric power transmission at the highest voltage 
with regard to different times of implementation of the corresponding 
measures (for example, with regard to temporary operation of new 
electric power transmission at reduced voltage): 3) selection of the 
directions and voltage of new electric power transmission, designed to 
generate the power of new power plants; 4) refining the location of nes 
power plants and substantiation of the plan of introduction of the 
generating capacities with regard to expenditures for development of the 
basic network of the electric power system; 5) substantiating the 
feasibility of introducing DC electric power transmission. 


It was necessary to solve a number of methodical problems during 
development of the program system. They were related to the need to 
increase the number of variables and of the number of calculating modes, 
to optimization of the flow distribution in a complex network, which 
combines DC and AC electric transmission, with regard to limitations on 
capital investments. Using the ORS E program system, one can simulate 
the process of development of energy associations in different 
Situations. Alternative measures on development of the network can be 
selected by using the program system, their efficiency for the entire 
range of anticipated future conditions can be determined with regard to 
space and time relationships, one can determine the factors that have a 
decisive influence on the optimal process of development, and one can 








determine which of the alternative measures and when they should be 
implemented to guarantee an optima! value of the optimization criterion. 
Integral criteria for a sufficiently long calculating period are used 
when compiling the strategic variants of development. This is related 
to the fact that the relative effectiveness for the more promising 
strategies during the initial period is lower, compared to extensive 
variants of development. The number of compared versions of development 
is so great that substantiation of the solution is impossible without 
special program systems for dynamic optimization of the development of 
electric networks. Due to the uncertainty of information, 
substantiation and decision-making are carried out on a dynamic plan: 

as the initial information is refined, the solution is also refined. 


The conditions of development includes information about the existing 
electric network, about the generating capacities and energy consumers, 
about future loads, about the calculating period, about alternative 
future measures, about restrictions, and about future technical and 
economic indicators. 


The concept of development step is used to describe the calculating 
period in optimization problems. This is related to the fact thet 
continuous progress of development of a system over time in practical 
dynamic calculation is always replaced by a discrete process. The time 
segment, during which the state of the system does not change, is 
considered a development step. The ORS E program system solves problems 
with number of steps up to 10. 


The input information for the ORS E dynamic program system is 
information about the calculating plan of the network and the calculated 
operating conditions of the network during the calculating period, and 
also about the set of possible measures on development of the network, 
representing the area of optimization itself. The spatial relations and 
area of search for the optimal variant according to development steps 
are simulated on the basis of the calculating plan of the network 
object. The main concepts that characterize the calculating plan are as 
follows: 


1) Notes. The combination of nodes of the network facility to be 
optimized is divided into groups (subsets), which are used to 
differentiate the technical and economic indicators. The nodes of this 
group are characterized by identical technical and economic indicators 
to compensate for the energy losses in the elements of energy 
associations. It is assumed in the ORS E when calculating the technical 
and economic indicators that energy losses in the elements of energy 
associations are covered in the outer nodes of these elements, i.e., in 
place. 


2) Elements of network. Consideration of the following types of 
elements of the network is permissible: a) of an AC electric 
transmission line; b) of a DC electric transmission Line; c) of a power 
plant (generators). The main notation of the element in the primary 








calculating scheme are the numbers of its outer nodes. Transformers are 
not directly shown in the calculating scheme. Lines with different 
rates of voltages may arrive at one node. The main technical and 
economic characteristic of the element of an energy association is the 
number of the set of technical and economic data. Standardization of 
elements is introduced by using the indicated number. However, 
additional economic indicators (constant components of capital 
investments) can be used, by which individual features of separate 
elements are taken into account. These indicators are also used to take 
into account the economic indicators of transformer and converter 
substations. 


3) Branches. The elements of the network are reflected in the 
calculating plan by branches. A branch is characterized by data of a 
specific element (AC or DC electric power transmission line, power 
plants) and by the number of parallel elements. Thus, a branch may also 
reflect several parallel LEP or generators with identical parameters. 
The primary scheme should contain both existing and future elements of 
the network. Thus, branches can be divided into existing and future 
branches. 


The number of elements according to development steps t can be 
increased. Parallel branches that reflect the elements of the network 
with different characteristics are permitted in the scheme. 


The connectedness of the calculating scheme is checked automatically by 
the programs of the system. The layout of the network is assumed 
connected if all its nodes, which can have a load, are connected to each 
other by existing or future branches. The characteristics which occur 
upon display of the electric network for dynamic optimization 
calculations include, for example, prohibition of elements with zero 
resistance (switches). One can require that one object be represented 
by several branches to display these measures. The calculating schemes 
of electric networks, compiled for dynamic optimization calculations, 
may have 100 nodes and 200 branches. Compilation of the calculating 
scheme for problems of optimizing the development of networks is a 
labor-intensive process and requires careful preparation of the input 
information even in the presence of check systems. But the labor 
expenditures to create the basic information about the calculating plan 
for optimization calculations of a given network are single- stage. 
This information is used repeatedly for a long period in studies of the 
development of EES. The calculating plan is only corrected during each 
study cycle. 


Compilation of the task for optimization (determination of the area of 
search for the optimal solution) and display of it in the calculating 
plan of the network predetermine to a significant degree the process of 
optimization. The erudition of the investigator is of important 
significance for compilation of the optimization task. A set of 
alternative measures, which are described by means of variables, is 
determined in compilation of the optimization task. The problems to be 











solved are those wi scre.@ variables. Onl, two values are possible 
for each variable: A measure corresponding to the variable is not 
implemented, and 1--the measure is implemented according to the 
variable. Each variable is characterized by an ordinal mumber and by a 
set of branches of the network, for which the number of elements of each 
branch per unit is increased on this and on subsequent steps of the 
optimization process after implementation of the variable at each 

step t. 


Wher postulating the optimization problem, to be solved by using the 
ORS E, an alternative measure should be determined such that 
implementation of it would yield an energy effect (reduction of losses 
in the network, reduction of loads and so on). 


The optimization program system for development of an ORS E uses 
variables of two types. Implementation of variables of the first type 
displays the measures on construction of electric power transmission 
lines over new routes, higher voltage lines and so on. No more than one 
new element appears in the corresponding future branches (with the 
number of elements in the primary scheme equal to zero). Implementation 
of variables of the second type shows the construction of new parallel 
LEP or installation of additional generators at electric power plants 
with the same parameters as existing plants, which is related to an 
increase of the number of elements both in existing and in future 
branches. When compiling the optimization task using the ORS E program 
system, one can impose restrictions on the possible period of 
implementation of variables, and the beginning and end of the period of 
implementation of alternative measures must be indicated. 


It becomes necessary in optimization problems to take into account the 
dependence between implementation of measures. For example, to replace 
the wires on some electric power transmission line, the line must first 
be constructed. Logical restrictions must be established for Boo} fan 
variables in these cases. There are two main types of thesog 
restrictions: 


{ 
‘ 


1) Dis,junction--states in which more than one measure (variable) from a 
given set is implemented are not permissible. Disjunction is used to 
simulate mutually exclusive measures on development of an electric 
network. This method of simulation can be used, for example, when 
selecting the optimal variant of the design version of an electric power 
transmission line (type of current, rated voltage, mark of wires). 


2) Limitations on the sequence of implementing measures to develop the 
network. If j after i is given, implementation of the j-th variable is 
not permissible if the i-th variable has not been implemented. 
Limitations of type "j after i” are used to simulate measures which can 
be implemented only in a specific sequence. For example, if temporary 
use of an electric power transmission line at lower voltage with 
subsequent conversion of it to rated voltage is considered, this 
conversion can be implemented only after implementation of the first 








measure--introduction of the electric power transmission line at the 
lower voltage compared to the rated voltage. 


The use of logical restrictions reduces the dimensionality of the 
dynamic problems. Thus, for example, the presence of 10 logical 
restrictions with 20 variables causes a reduction of the number of 
possible states by approximately a factor of seven. Compilation of the 
optimization task is made difficult by the circumstance that the number 
of variables in dynamic optimization models is limited: the maximum 
permissible number is 64 variables. Accordingly, some aggregation 
(consolidation) of measures is possible. 


Let us consider below the methods of simulation and optimization which 
are used in the ORS E program system. 


Notations 

q --number of nodes in calculating plan; 

n --total number of branches in calculating plan; 

non --total number of branches that show electric power 
transmission line (LEP) ; 

Ngoc --total number of branches that show electric power 
plants (ES); 

i --ordinal number; 

c(i, t) -~-number of parallel elements of i-th branch on t-th 
development step, 

for LEP--number of parallel electric power 
transmission lines, 
for ES--number of units of same type; 

r(t) --number of calculated sets of loads, considered on 
t-th development step 

1(i) --length of electric power transmission line for i-th 
branch, km; 

Pi, t, J) --flow of active power in i-th branch on step t for 
i-th calculated load, MW; 

R(i) --coefficient proprotional to effective resistance, 
for determination of active power losses in one 
element of i-th branch, 1/MW; 

X(i) --coefficient proportional to reactance, for 
determination of reactive power losses, 1/MVA; 

Ygfi) --coefficient proportional to conductivity, for 
determination of flow of active power in one 
element of i-th branch (economical conductivity of 
branch 1), MW; 

Xd --longitudinal reactance of generator, percent: 

PL (i) --specific active power losses to corona in i-th 


branch, kW/km; 
cos --power coefficient of power flux, relative units; 


of 








g.(1) --specific susceptance of electric power Lransmission 


P | 
line, 10 ° com/km; 

Uti) --rated voltage in i-th branch, k\; 

K i) --capital investments in LEP, shown by i-th branch, 
dependent on its length, thousand rubles/km; 

Keg !t) --capital investments in i-th branch, independent of 
its length (including those for substations), 
thousand rubles; 

NOCT |. | _ a, . 

k gc!) --capital investments in i-th branch, shown by 
electric power plant, independent of number of 
units, thousand rubles; 

Key (2) --capital investments for installation of one unit in 
i-th branch, representing power plant, thousand 
rubles; 

Pen --capacity of unit, MW; 

Coonan --specific fuel expenditures, kopecks/kW-hr; 

Take | i) --number of hours of utilization of installed power of 
i-th power plant; 

Ky --specific capital investments for installation of 
reactors, thousand rubles/MVA; 

K. --specific capital investments in pressurizer on one 
electric power transmission line, thousand 
rubles/MVA; 

”» --degree of compensation of charging power of power 
transmission lines, relative units; 

' --efficiency of capacitive power of electric 
transmission line, relative units; 

Ava --relative costs, dependent on capital investments 
Keg!) for electric power transmission lines, 
relative units; 

a0 --relative costs, dependent on capital investments 
Kea (1) for power plant, relative units; 

a, --relative costs for depreciation and maintenance for 
electric power transmission lines, dependent on 
K, | i), relative units; 

A OCT --relative costs, dependent on capital investments 


Ke gal i), relative units; 











Kpeak --coefficient for determination of expenditures for 
compensation of load reactive power losses, thousand 
rubles/year :-MVA; 


340) t) --relative expenditures for losses in i-th branch, 
dependent on loads, kopecks/kW-hr; 

Ree t) -~-relative expenditures for no-load losses in i-th 
branch, kopecks/kW-hr; 

Tr(J) --length of j-th mode, hr; 

3, --dynamic reduced expenditures; 

346 --additional criterion for estimation of disruption 
of power balance conditions; 

3, --additional criterion for estimation of violation of 
restrictions on power flows; 

kk --number of parts of network unrelated to each other: 

Pui, J) --power unbalance in j-th part of network for i-th 
set of calculated loads, MW; 

Rt --hypothetical resistance; 

L(j) --set of branches representing LEP for which 
restriction on power fluxes is disrupted, at j-th 
set of calculated loads; 

Pron! ?? --capacity per circuit, MW; 

A and k --coefficients for calculation of additional criterion 


for estimation of violation of restriction on power 
fluxes are given in input information; 


6 --number of years during which capital investments are 
made or annual costs vary; 

Ey --normative factor of comparative effectiveness of 
capital investments; 

Eun --normative factor of reduction of expenditures at 
different times; 

3.10) --reduced expenditures throughout system during vear 
v; 

3.19) --reduced expenditures throughout system at normal 
operating step; 

v, --year of reduction of expenditures at different times 


(the first vear of the calculating period is used as 
the year of reduction in the ORS E program system); 


I(t) --first year of t-th phase; 

S(t) --coefficient of specific weight of expenditures 
3. (dtd) on t-th development step in specific 
function; 

e(t) --state of network on step t (combination of 


implemented measures on development of electric 
power system during time up to t, inclusively); 
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fit, e(t)) --minimum value of functional--specific function, 


corresponding to optimal process of development 


during t steps from initial state to some state 

e(t); 

git, e(t)) --component of functional, dependent only on step t 
and state e(t), but independent of path of 
transition to given state e(t); 

e(t, m) --states of development on step t with m implemented 
measures ; 

f 


--constant, called critical gradient value; 
grad*f(t, e(t, m))--relative maximum increase of functional for state 


e(t, m) according to comparison to functional for 
state e(t, m- 1) (relative gradient). 


Indicators and Criteria of Variants 


Active power losses AP(t, j) on step t with the j-th set of calculated 
loads are determined as the sum of load losses and no-load losses: 


AP (t, |) =APul(t, j) +APxx (t, j). (1) 


The active power load losses in MW are 








cae: eel 
sPu(t.i)= Dy — | (2) 
oss c l,i 


The active no-load power losses in MW are: 


non 
~ 
APxx(t, j) = » Py(i) -10-3-¢ (i, t) «1 (i). = 


The reactive power losses in the electric power transmission lines in 
MVA are : 





~~ P;- it. pyX) 
\Q(t, j)= » ( 7h ine ae 


— 5  Bvli) Ua? (i) 10-8 e(3, £) -1(i) ) (4) 
~ cos’ qg-cté, ! 


The annual costs are determined as the sum of the costs for depreciation 


and maintenance, for energy losses and for compensation of the reactive 
power in thousand rubles/year: 











H(t) = Wi) + a(t) +H (ty. 


The annual costs for depreciation and maintenance in thousand 
rubles/year are: 


"non 


(t) = py Hi: nan (i, t) + y Mx. ac (i, ft). 


The annual costs for depreciation and maintenance for branch i, 
representing the LEP, in thousand rubles/year are: 


Mx, nan (i, t) = Kya (i) ( (C(t, f) —C(i, 0) )anat+ (Ka(i)-aa+ 
+ Zp(t) - Uy? ( (¢)-10-°(K, ‘Yp— Keyc) axa) l (i) ( (c(i, t) —c(i, 0)). 


The annual costs for depreciation and maintenance for branch i, 
representing the LEP, in thousand rubles/year, are: 


Ux, 9c (i, t) =Koali)- (c(i, t) —¢ (i, 0)) -aoat 


gDOCcT 


+(Kgc (i) + Kxa (i) ) - Qnoer. 


The annual costs for losses in thousand rubles/year are: 


P2(i,t RW. in. torte 
H,(t)= y Cli. t) °3 L(t, t) T-({) 10 : 








+ ¥ Py (i) -3/7 (i, t) <1 (i) -8760- 10-* (c(i, £) —c(i,0)), 


= 


The annual costs for compensation of reactive power in thousand 
rubles/year are: 


"n9n 


K+ GyaX (i 
Hg (t) = py — max P,? (i, ¢, j). 





cos? — 


(6) 


(8) 


(9) 





‘apital investments are determined as the sum of capital investments in 
the electric power transmission line and the power plant: 


"nan 


"3c 
K:(t)= y Kaan (i, t) + » Koc (i,t). (10) 





Capital investments for branch i, representing the LEP, in thousand 
rubles, are: 


Kasn (i, £) = (Kua (i) + (Ka(i) + gp (i) -Un?(i) - 10> 
. (Kpyp — Key, ))-L(i)}- (c(t, f) —c(i, 0)). (11) 


Capital investments for branch i, representing the power plant, in 
thousand rubles, are: 


K3-(i, t)= Kec +Kea(i) -(c(i, t) —c(i, 0)). (12) 


Dynamic expenditures are determined by the expression in thousand 
rubles/year: 


> 
33> )) (Ea-Ke(t) + H(t) -S (2). (13) 


The dynamic reduced expenditures for the step process (the calculating 
period is divided into t normal operating steps, i.e., the loads for all 
calendar years of the step t are equal, while the scheme and parameters 
of the system do not vary during the step) are determined by the 
expression: 


4 
3.= ) 3. ((t)) S(t) +3. (0 - 1) 1 +Enn)*%e-?, (14) 


where the coefficient of the specific weight of expenditures on the t-th 
development step is determined by the formula: 


t 


Si: »)® Are (1) + Ew) °e-*'* (15) 











The ORS E program system permits the investigator to assign the value of 
coefficients S(t) for t = 1, 2, ..., T according to his discretion. If 
the values of S(t) are not assigned, they are determined by formula 
(15). An additional criterion of estimating violations of the condition 
of power balance on step t are determined by the expression in thousand 
rubles/year: 


qt 


KR 
Zac (t) = y Py? (i,j) -R:-Aca:T, (i). (16) 


The additional criterion of estimating violation of restrictions for 
power flows on step t is determined by the expression in thousand 
rubles/year: 


rt) I 
Z(t) = ¥ »y (Ex-Kaan (i. t) — Vy aon. t))> 


p ne? (17) 
— absPrii,t, i) \) 


—— 





rit) (i, f)- Pant) 


Model of Power Flows 
The flow distribution of active power in the ORS E program system is 
determined by the minimum power losses in the network with regard to 


Kirchhoff’s first law as supplementary conditions |1]. The expression 
for optimal power flows has the following form: 


P,(l) Vf) (8) &,( )), (18) 


where b,(i) and bt }) are the economic potentials for the initial and 


final node of branch |, respectively. The economic potentials are 
determined by a system of q equations: 


Y,-A,=P. (19) 


where A, is the vector column of the nodal economic potentials, P is the 


vector column of the loads of the nodes, and Y, is the quadratic matrix 


of economic conductivity. 





The active power losses in transformers in system-forming networks are 
relatively low compared to losses in electric power transmission lines. 
Therefore, the calculated layouts of the network do not contain 
transformer branches when determining the optimum flow distribution (see 
figure). Caiculation of the optimal flow distribution permits 
convenient simulation of parallel operation of an AC network with DC 
elements [1]. 
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Calculated Layout of Optimization of Development of 
Energy Association: 
l--existing electric power transmission lines (LEP) of 500 kV; 
2--existing LEP of 500 kV, temporarily using 220 kV; 3--existing 
branches with several elements; 4--future and alternative LEP of 
500 kV; 5--future and alternative LEP of 1,150 kV; 6--future and 
alternative DC LEP of 1,600 kV; 7--future and alternative power 
plants 


KEY: 
1. Regional power plant 2. Existing elements 


Target Function 


The following target function is used when searching for the optimal 
variant of development: 








J Fa3at Zu +Zo. (20) 





The supplementary criterion of estimating the violations during the 
entire calculating period is determined by the expression: 


\ (21) 
/ Z y 2()S() 


The ORS F program system txukes into account restrictions on capital 
investments Kft) according to development steps t = 1, 2, ..., T. If 


the values of K A(t) are mot given, it is assumed that K,{(t) - 109 rubles 


for : Be eseg Oe 
Me .1 pt imi zation 
All eps t = 1, 2, ..., T are considered sequentially during 


yptimi zation. ‘The optimal development during t steps to state e(t) is 
jetermined at each optimization step [2]. To do this, the following 
recurs on relation is used: 


j(t,e(t))=g(t, e(G) + min f(t—1,e(t—i)}. 
(t=-1)er(t) (22) 


Oniy ptimal initial states, which are local optimums according to the 
following condition, are sought during the next optimization step |t 


P(t, ety) min J(f,e (f)). (23) 


e§ 


The optimal initial states (OIS) and the corresponding values f(t, e(t)) 
are stored, which also permits one to determine the optimal development 
on step t + 1 during t + 1 steps to any state e(t + 1). The optimal 
initial states on the next step t is sought by sequential consideration 
of the subset of states {e(t, m)}, m= 1, 2, ..., where m is the number 
of implemented measures. 


Because of the large dimensions of the problems of optimizing the ORS E 
program system, the approximate method of finding the optimal initial 
state is used. 











In this case, the optimal initial states are selected by the feature 


—grad*/(t, e(t,mj)) <6. (24) 


The relative gradient is determined by the formula: 


f(t,e(t,m))-— min f(t, e(t,m—1)) 


ett, acm liee(t, m) 
grad f(t, e(t, m))= — A -100%. a 
[(t,e(t, m)) 








The search for the optimal initial states is stopped if no optimal 
initial states are found at given m. One could find all the local 
minimums of the target relation at é = 0 and upon retention of all 
optimal initial states at each optimization step and one could find the 
global optimum. However, this is unrealistic due to the large 
dimensions of the dynamic problems. The real dynamic optimization 
problems can be solved if the number of optimal initial states (6 > 0) 
is retained at each development step by retaining only the optimal 
initial states with the best value of the functional for different 
competitive strategies. To reduce the set of optimal initial states in 
the program, one can also use the restriction form above FAA(t) for each 
step t = 1, 2, ..., T. Only those optimal initial states e’(t) for 
which the following condition is fulfilled are retained during 
optimization 


Ie) erg aity. (26) 


min f(t, e(t)) 





If FAA(t) is not given, it is assumed that FAA(t) = 1.4, since 
experimental and theoretical studies show that, if the values of the 
functional of some state differ from the minimum value more than 1.4, an 
optimal strategy may not pass through this state. 


Pseudo-dynamic optimization is realized at FAA(t) = 1. Thus, the values 
of FAA(t) can be given in the range 1-1.4. The calculating time 
decreases as the value of FAA(t) decreases, but the probability of a 
loss of optimal strategy increases. The values of FAA(t) for subsequent 
optimization steps can be assigned less than for the initial values. 

The maximum possible number of retained optimal states is WM = 40. A 
decrease of WM reduces the time for optimization, but the probability of 
a loss of the optimal variant according to the integral target relation 
during the calculating period increases. Pseudo-dynamic optimization is 
realized at WM = 1. The procedure of optimization (22) for steps t = l, 
2, ..., T is performed repeatedly (by the iteration method) at 4 > 0. 
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Table 1. 


Optimized Measures (Variables) 

















| (1) ( 2) MHOMECTBO BeTBCII | (3 Prpakusenna 

i | HuuMenosanue MeponpuaTHa | (L}. | (L), | (LO}, | (tq—ty), | 

| | 2 | * H §° | b* 

1 Bsor s2ekTpoctaHuNH B vate 14 {4} 31 — -- ‘2-2 - 

2 Baor sx1eKktpoctanuiun a v3te 15 5) 32 — -- 2—2 — 

3 Bsor s2eKkTpoctasuNkK B V3te 3 6 33 -- _- 2—?2 

4 Baoa 13M «xB 12—14 7) 20 — — I—3 — 

5 Bsoa .1911 500 KB 12—15 8) 21 — — J—3 -- 

6 Bsor asvxuennoit .139T1 500 KB 12—13(9) @ — — I—3 — 

7 Beoa 11911 500 KB 3—12 (10) 23 . — -—- I-33 = 

8 Bsoa .1311 500 KB 10—11—12 (33) 24,26 — — l—J — 

9 Bsoa -1311 500 kB 6—11—12 (12) 24,25 — -- J—3 — 

10 Bsoz .13M 500 «KB 13—14 (13) 19 -— 23 

11 Bsoa 13M 500 xB 13—15 ( 14) 18 — — 2—3 — 

12 Bsoa 3-ii napasaeabnoi £1311 500 KB 1-16(15}+ 33 — i—3 ~ 

13. Bao 4-ii napaaseapuoit £1311 500 kB 1--1€ 16+ = 3 — 1—3 12 

14 Bsog 5-4 napaanerabxoli .1311 500 kB I-I@ 1A OSC 1—3 13 

15 Bsoa 3-ii napaaszeavnolt £131 500 xB (18) —- 9%6 — i—3 - 
l—9—3 ; 

16 Bsog 4-i napasseswxol .1311 500 KB (14) - 96 — 1—3 15 
3 _ 

17 Besos 5-it napaaneapnol .13M1 500 KB (20) - 3$6&6 - iI—3 16 
1—2—3 

18 Besos 2-ii napasreanoii .1311 500 xB (21) — 910 - | ang} — 
3—5—6 

19 Seon 13911 1150 KB 6—16 (22) 28 -- _ I—3 — 

20 Baoa .13911 noctrosxunoro toxa 6—16 (333 9 —_ -- I—3 = 

2! Tlepesoa JISII ¢ 220 na 500 KB (2 297 _ 17 I—3 — 

Notes: i--number of variable; {L},'--set of future 
branches for which the number of elements assumes the 
value of 1 after implementation of the i-th measure; {L},°-- 


set of existing and future branches for which the number 
of elements increases by 1 after implementation of the 


i-th measures; 


{Lo}, --set of branches for which conductivity 


is set equal to zero (the branch is disconnected) after 
implementation of the i-th measure; 


permissible implementation of i-th measure (t 


number of the optimization step); i 
after implementation of which implementation of measure i 


is permitted. 


Name of measures 

Set of branches 
Restrictions 

Introduction of power plant 
Introduction of power plant 
Introduction of power plant 


(Key continued on following page) 


in node 
in node 
in node 


3¢ 


(t ),--period of 


—t 
H K 
is the ordinal] 
--number of measure, 














7. Introduction of LEP 12-14 of 500 kV 

8. Introduction of LEP 12-13 of 500 kV 

9. Introduction of two-circuit LEP of 12-13 of 500 kV 
10. Introduction of LEP 3-12 of 500 kV 

11. Introduction of LEP 10-11-12 of 500 kV 

12. Introduction of LEP 6-11-12 of 500 kV 

13. Introduction of LEP 13-14 of 500 kV 

14. Introduction of LEP 13-15 of 500 kV 

15. Introduction of third parallel] LEP 1-16 of 500 kV 
16. Introduction of fourth parallel LEP 1-16 of 500 kV 
17. Introduction of fifth parallel LEP 1-16 of 500 kV 
18. Introduction of third parallel LEP 1-2-3 of 500 kV 
19. Introduction of fourth parallel LEP 1-2-3 of 500 kV 
20. Introduction of fifth parallel LEP 1-2-3 of 500 kV 
21. Introduction of second parallel LEP 3-5-6 of 500 kV 
22. Introduction of LEP 6-16 of 1,150 kV 

23. Introduction of DC LEP 6-16 
24. Switch of LEP from 220 to 500 KV 


The optimization is performed on the first iteration at some given value 


Cay and considers all variables. The value of é decreases on 


subsequent iterations at given step Aé (Aé < Oey: Preoptimization of 


the optimal variant of development is performed. fhe ORS E program 
system presents the investigator the opportunity to perform 
preoptimization of some other variant as well if desired. In this case, 
the investigator should indicate the ordinal number of the preoptimized 
variant. If the ordinal number is not given, the optimal variant is 
preopt imized. 


Thus, the optimization problem is solved with a reduced number of 
variables and with increased accuracy of finding optimal initial states. 
Optimization ends when the values of 4 reach zero, i.e., when finding 
the exact values of the local minimums of the target relation, 
corresponding to the sought optimal initial states, it is guaranteed 
during searching for the optimal initial states. 


Example 


The calculating scheme for optimization of development of a system- 
forming network with refinement of the location of new power plants is 
presented in the figure. The considered energy association is separated 
into two regions--RES 1 and RES 2, having different cost of energy. The 
period to be studied is divided into three development steps with length 
of 5 vears. A tota! of 10 calculated sets of loads is considered on 
each development step. Introduction of power plant 21 is provided 
during the first optimization step, and the problem of generation of 
power by this station must be solved. Another power plant must be 
introduced on the second step, and its location must be refined (node 
14, 15 or 3). The main problem to be studied is that of selecting the 
variant of electric power transmission from node 6 to node 16. There 








are three alternatives: 1) enlargement of the 500-kV network, 2) 
construction of a 1,150-kV LEP and 3) construction of a DC LEP. A 
special problem which is also studied in this problem is determination 


of the moment of transition of LEP 8-17, temporarily utilized on 220 and 


5900 kV (construction of a 500/220-kV substation at node 8). 


Table 2. Logic Restrictions of Type “Or” 


Set of variables from which 
implementation of no more than 
one variable is permissible in 

Number of restriction one development variant 








1, 2, 3 


whore 
U 
Oo 


Table 3. Results of Optimization 

















(1) (2 )CPOKH peaH3aulHH 
Homep nepemeHHoA i 7 | _ 
| 2° | 3° 
l l l 
9 — — 
3 2 2 
4 1 
5 ] ] 
6 — iam 
7 — _ 
8 l l 
9 '] — 
10 — — 
1] — ~— 
12 - 3 
13 — — 
14 — —_ 
15 - ome 
16 - - 
17 - — 
18 2 2 
19 3 3 
0 — — 
21 l l 
*Notes: t,!--ordinal number of optimization step 


at which measure i was implemented in optimal variant of 
development; tj ’?--ordinal number of optimization step in 
variant for which the value of the functional is 1 percent 
higher than the minimum value. 


1. Number of variable i 2. Periods of implementation 
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A list of measures to be optimized according to development of the 
network is presented in Table 1. Variables 12-18 are those of the 
second type. They are described by a set of existing branches, for 
which the number of elements is increased by one by implementation of 
the measure. Variables 15-18 are described by two branches, while the 
remaining variables of second type are described by only one branch. 

The specific sequence of realizing them is observed when using variables 
of the second type and logic restrictions ‘after” are used (Table 1, 
column 7). Restrictions for the period of implementation of alternative 
measures (Table 1, column 6) and logic restrictions of the "or" type 
(Table 2) are also used in the considered example. 


Introduction of these restrictions is determined by postulation of the 
optimization problem. For example, those layouts of the network are 
impossible in which branches 21 (variable 5) and 22 (variable 6) are 
present simultaneously, since these branches represent the same line in 
a one- or two-circuit version. Implementation of more than one measure 
from set 1, 2, 3 is also impossible, since only one new power plant is 
planned for construction during the second development step. 


The results of optimization are presented in Table 3. 
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